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^ I Gravitational microlensing provides a new technique for studying the 
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c^ I surfaces of distant stars. Microlensing events are detected in real time and can 
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be followed up with precision photometry and spectroscopy. This method is 
particularly adequate for studying red giants in the Galactic bulge. Recently 
we developed an efficient method capable of computing the lensing effect 



*> I for thousands of frequencies in a high-resolution stellar spectrum. Here we 

j^^ I demonstrate the effects of microlensing on synthesized optical spectra of red 

^ I giant model atmospheres. We show that different properties of the stellar 

O I surface can be recovered from time-dependent photometry and spectroscopy of 



a point-mass microlensing event with a small impact parameter. In this study 
we concentrate on center-to-limb variation of spectral features. Measuring such 

Ph| variations can reveal the depth structure of the atmosphere of the source star. 

O 

> 

^ ' Subject headings: gravitational lensing — red giants — stars: atmospheres 

Submitted to The Astrophysical Journal, February 1999 



-2- 



1. Introduction 



Stars populating tlie red giant branch spend a substantial fraction (> 5%) of their 
entire hfetime there. Red giants are thus a major component of any galactic stellar 
population. In old systems, like Galactic globular clusters, their contribution is even 
larger. The parameters of the red giant branch have been traditionally derived from 
multi-color photometry and interpreted using stellar atmosphere models (e.g. Vandenberg 
& Bell 1985; Kurucz 1992). These model atmospheres involve one- dimensional integrations 
of semi-infinite slabs in hydrostatic and radiative equilibrium. However, real atmospheres 
of red giants are usually extended and the slab approximation is not adequate. Despite 
significant progress in computing extended spherically-symmetric models (e.g., Scholz & 
Tsuji 1984; Plez, Brett, & Nordlund 1992; to name a few), it is still difficult to use these 
improvements in massive grids of models, which therefore continue to be computed as 
plane-parallel (Houdashelt et al. 2000). A number of other assumptions and simplifications 
had to be usually made, some of which have not been verified yet by direct observations. 
Examples include the treatment of convective transport and small-scale velocity fields (e.g. 
micro- and macro-turbulence). Similarly, stellar evolution models also depend on a number 
of parameters which describe poorly understood physics, such as convection near the surface 
of the star. 

These problems have hindered significant improvement in the determination of 
fundamental stellar properties over the past two decades. The angular diameters of only 
a small number of stars have been measured up to now, with an accuracy of ~ 10% at 
best (Armstrong et al. 1995). Stellar disk brightness distributions (limb darkening) have 
been inferred only for a handful of stars at the level of a proof-of-concept (directly — 
Mozurkewich et al. 1991; or in binaries — Andersen 1991). An exception is provided by the 
direct HST observation of the red supergiant a Ori, which revealed limb darkening as well 
as a bright region on the surface (Uitenbroek, Dupree, & Gilliland 1998). In general, the 
available inferred results cannot be used to build a model atmosphere, as in the case of the 
Sun. Due to this lack of an adequate observational basis, our understanding of stellar light 
is implicitly tied to the solar atmosphere model. 

This situation is particularly unfortunate, since recently there has been an increased 
demand for accurate stellar models in a number of applications. For example, much 
improved color-temperature conversions are needed for calibrating distance indicators and 
determining stellar ages. Stellar population syntheses (e.g. for high-redshift galaxies) also 
rely critically on the accuracy of current stellar models. Our Sun is not a good standard for 
most of these applications. 
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The advent of large interferometric arrays (Keck interferometer, CHARA, SIM, VLTI) 
will go a long way towards solving these problems, but these complex and expensive 
facilities are still mostly in a preliminary development phase. In the meantime, gravitational 
microlensing offers an easily accessible, immediate, and inexpensive means for imaging at 
least large stars, such as red giants. It also offers, by its nature, access to stellar populations 
in the Galactic bulge and the Magellanic Clouds, which are well beyond the reach of any 
interferometer. 

By now there are more than 350 microlensing events detected towards the Galactic 
bulge and the Large and Small Magellanic Clouds (see e.g. Alcock et al. 1997a, 1997b, 
1997c; Becker et al. 1998; Renault et al. 1997; Afonso et al. 1999; Udalski et al. 1997). 
During a Galactic gravitational microlensing event the flux from a distant star is temporarily 
amplified by the gravitational field of a massive dim object passing in the foreground. The 
standard light curve in the point-source limit is characterized by two observables: its peak 
amplitude and duration (Paczyhski 1996). These observables physically depend only on the 
lensing impact parameter and the Einstein radius crossing time of the source. However, 
in events with a small impact parameter the light curve is also affected by the resolved 
surface of the lensed star, thus providing an excellent opportunity for stellar surface imaging 
(Valls-Gabaud 1995, Sasselov 1996). 

The first well documented case of resolved finite-size microlensing effects is MACHO 
Alert 95-30 (hereafter M95-30), with the impressive follow-up campaign by the GMAN 
and MACHO collaborations (Alcock et al. 1997d). The first spectroscopic observations of 
a binary lens event in which the caustic crossed the face of the source star were reported 
by Lennon et al. (1996). In a similar caustic-crossing binary event, Albrow et al. (1999) 
recently determined limb darkening profiles of the lensed K giant star for both observed 
spectral bands (/ and V). 

Finite-size effects were first theoretically studied as methods to partially remove the 
degeneracy of microlensing light curves. These methods are based on alterations of the 
standard light curve (e.g. Gould 1994, 1995; Nemiroff & Wickramasinghe 1994; Witt & 
Mao 1994; Gould & Welch 1996), resolved polarization (Simmons, Willis, & Newsam 1995; 
Simmons, Newsam, & Willis 1995; Bogdanov, Cherepashchuk, & Sazhin 1996), spectral 
shifts due to stellar rotation (Maoz & Gould 1994), and narrow-band photometry of 
resonance lines (Loeb & Sasselov 1995). Here we want to put emphasis on the inverse 
problem — using microlensing for studying stellar surface features and probing the 
atmosphere of the source star. 

The resolved stellar surface brightness distribution, B{X,f), can vary strongly with 
wavelength in selected spectral regions - in the continuum as well as within spectral 
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lines. The time-dependent microlensing amplification then becomes wavelength-dependent 
through B{X,f). Photometry in different bands and sets of spectra taken in the course 
of a microlensing event can therefore be used to study the brightness distribution of 
the source. Microlensing light curve chromaticity was first evaluated by Witt (1995), 
Valls-Gabaud (1995) and Bogdanov & Cherepashchuk (1995). Methods for obtaining 
the brightness distribution by light curve inversion were studied by Bogdanov & 
Cherepashchuk (1996) and Hendry et al. (1998), inversion error analysis was presented 
by Gaudi & Gould (1999). Few works have been published so far on spectral effects in 
stellar microlensing. These include studies of line profile changes in rotating giants (Maoz 
& Gould 1994; Gould 1997) and in circumstellar envelopes in bulk motion (Ignace & 
Hendry 1999). First results for selected spectral lines in a simple model stellar atmosphere 
were presented by Valls-Gabaud (1996, 1998). 

In order to optimize observations of spectral effects in actual events, the most practical 
approach is to theoretically predict the lensing effect on the entire synthesized optical 
spectrum of the source in microlensing alerts similar to M95-30. This way one can determine 
the most sensitive spectral features before the event peaks. An efficient method and 
computer code capable of dealing with the tens of thousands of frequencies involved was 
developed by Heyrovsky & Loeb (1997; hereafter Paper I). In this work we use the method 
to explore point-mass microlensing effects on both low and high resolution optical spectra, 
synthesized directly from up-to-date red giant model atmosphere calculations. 

The outline of this paper is as follows. In §2 we describe the model atmospheres used 
in this paper. The effect of microlensing on the overall spectrum of the source star as 
well as on individual spectral features is presented in §3. Finally §4 summarizes the main 
conclusions and future prospects of this work. 



2. Model Atmospheres of Red Giants 

In a stellar population with some spread of metallicities red giants exhibit a range of 
temperatures. The cooler giants have atmospheres which, unlike in the case of the Sun, are 
dominated by molecular compounds. Most of the hydrogen gets locked into H2 and most 
of the carbon into CO. As the temperature drops, TiO, VO, and eventually H2O become 
prominent. This requires paying special attention to the equation of state and background 
opacities (Hauschildt et al. 1997; Kurucz 1992). Of great importance is the treatment of 
convection, as transport by convection affects a large part of the photosphere. Finally, for 
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Fig. 1. — NLTE line profiles from different positions on the stellar disk in a red giant model 
with T=3750 K, log5'=0.5 . Left panel: profiles of the hydrogen Ha line from disk center 
at /i=0.996 (lowest curve), through /i=0.5 (central curve at line core), to very near the limb 
at /i=0.004 (upper curve at line core). Here /x is the cosine of the angle between the normal 
to the stellar surface and the line of sight (referred to as "the incidence angle" hereafter). 
Right panel: same for the hydrogen H/5 line; the profile sequence is /i=0.996 (lowest curve), 
/i=0.004 (central curve) and /i=0.5 (upper curve). Note the difference in profile shape, 
as well as in its change with /i. 



models with high luminosity (and low surface gravity, log g < 1), spherical geometry must 
be adopted in the treatment of radiative transfer (Hauschildt et al. 1999). 

Here we have limited our study to static (i.e., non-Mira) red giant models. Our main 
interest is to compute realistic models of the center-to-limb variation of the intensity on 
the stellar disk — both for individual spectral features in high resolution, as well as for the 
optical pseudo-continuum. The state-of-the-art work for such models has been done recently 
by Hofmann & Scholz (1998) and Jacob et al. (1999). They employ model atmospheres 
from the grid of Bessell et al. (1989, 1991). Here we use similar models which reproduce 
their results. Striving for this level of sophistication is crucial in our study of center-to-limb 
and depth dependence of spectral features in red giants. For example, Valls-Gabaud (1998) 
used the linear limb darkening law and fitted simple analytical expressions for spectral 
lines from Kurucz (1992) model grids to predict chromatic and spectroscopic signatures of 
microlensing events of stars in general. We could not use the same simple approach, because 
we find that for red giants these linear approximations fail by a factor which is much larger 
than the uncertainties of the models, in line with similar findings by Jacob et al. (1999). 

The structure of the models in this paper is derived by assuming local thermodynamic 
equilibrium (LTE) for the background opacities in flux-constant plane-parallel atmospheres. 
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Fig. 2. — Line-forming regions of the hydrogen Ha and H/5 hnes shown in Figure |I|, 
represented by their contribution functions (CFs). The CF contours range from 0.1 
(outermost) to 1.5 (innermost) in units of normahzed intensity. Left panel: the hydrogen Ha 
hne; right panel: the hydrogen H/5 line. Zero height is defined at continuum optical depth 
of unity for the rest wavelength of each line respectively. 



We use the most recent opacity sampling routines of Kurucz (1999) and TiO line data of 
Schwenke (1998). These models serve as input for computations of hydrogen and calcium, 
which are treated out of LTE (NLTE) with multi-level model atoms and in spherical 
geometry (for more details see Loeb & Sasselov 1995). The hydrogen atom model has 
energy levels n=l— 5, plus continuum; all bound-bound and bound- free transitions are 
calculated in detail. The calcium atom model contains 8 energy levels for the neutral species 
(Cal) and 5 energy levels for the first ionization state (Call), plus continuum (Calll); all 
transitions are calculated in detail. In the case of red giants, atomic hydrogen and calcium 
are not important for the overall atmospheric structure. However, they still provide strong 
spectral features with significant center-to-limb variation (e.g. in the hydrogen Balmer lines. 
Call H&K, etc.). Most moderately strong and weak lines diminish towards the limb (see 
Figure |l|). As discussed further below, this variation could be detected in a microlensing 
event by measuring the total equivalent widths of spectral lines on medium-resolution 
spectra with a high signal-to-noise ratio, or it could be seen directly on high-resolution 
spectra (see §3). 

The variation of the intensity from the center of the stellar disk to its limb is studied 
theoretically by computing the contribution function of each transition. This function 
defines the line-forming region in both space (depth in the atmosphere) and frequency. 
The line-forming region is shaped by local as well as non-local influences. Following 
Magain (1986), we define the NLTE contribution function (CF) for the relative line 
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depression as 

where /i is the cosine of the incidence angle; r^p is the optical depth at a reference wavelength 
Ao ; tr is the optical depth corresponding to the opacity hiR = ki + Kc Sc/h ', f^o is the 
absorption coefficient at Aq ; ki and Kc are the line and continuum absorption coefficients, 
respectively; Si and Sc are the line and continuum source functions, respectively; and Ic is 
the emergent continuum intensity (if the line were absent). The traditional dependence on 
frequency has been omitted from our notation for clarity. The CF for a given synthesized 
spectral line indicates where the line is formed, and thus it provides a theoretical tool for 
recovering the depth dependence of thermodynamic variables in stellar atmospheres. 

The CFs of the two lines presented in Figure ^ are shown in Figure 0. Note how two 
spectral lines from the same element and spectral series can exhibit differences in their 
CFs, in the CF's location against the local continuum, and hence in their center-to-limb 
variation. For example, the change in the strength of the H/? line undergoes a reversal 
close to the limb of the star. The line gradually weakens away from the disk center until it 
reaches a minimum, then it becomes slightly more prominent again near the limb. In the 
case of the Ha line this effect is much weaker and occurs closer to the limb. This difference 
between the two lines is due to the stronger dependence of H/3 on the TiO opacity, which is 
dominant at its wavelength. 



3. Effects of Microlensing on the Source Spectrum 

3.1. Microlensing of an Extended Source 

Next we study the effect of microlensing by a point-mass on the spectrum of the source 
star. For simplicity we describe the lensing configuration in terms of angular distances in 
the plane of the sky, using the angular source radius as a length unit. 

In the absence of a lens, the observed flux from an unresolved source star at 
wavelength A is obtained by integrating the surface brightness distribution -B(A, r ) over the 
projected surface of the star S, 



F^{\) = j B{\r)dT. . [1] 



Any information about the surface structure is therefore concealed by the integraL 
This information can be uncovered by microlensing, as its signature contains differential 
information on the brightness distribution. The flux from a simple point-source separated 
from a point-mass lens by a projected distance a in the plane of the sky gets amplified by a 

factor 

, , , a'^ + 2e^ 

where e is the angular Einstein radius of the lens (Paczyhski 1996). The lensed flux from 
an extended source can then be computed as 

F{X,ao)= jB{X,f)Ao{a)dJ: , (3) 

s 

here a is the distance from the lens to point r on the source, (Tq is the distance to the 
source center. In our units, (Jq = 1 corresponds to having the lens projected at the source 
limb. The conversion between the a distances and the source-centered position vector 

r = {p cos ip , p simp) is 



a = ycTQ - 2aop cos ip + p"^ . (4) 

The coordinates are oriented so that the lens lies in the direction ip = 0. In this notation 
p = \/l — /i^, where p is the cosine of the incidence angle used in §2. Generally the flux in 
equation (^ depends also on the angle between the lens and a fixed direction on the source 
(see Paper I). However, in this work we consider for simplicity a circularly symmetric source 
with B{\,r) = B{X,p). Angular surface irregularities such as spots are discussed in detail 
in a separate paper (Heyrovsky & Sasselov 2000a). Equation (|]) indicates that as the lens 
moves, it gradually scans the source, giving highest weight to the region of the star closest 
to the projected lens position. The total resulting amplification of the source star is 

A{X,ao) = ^r^ • (5) 

The center-to-limb surface brightness profile {B as a function of p) at different 
wavelengths has not only a different amplitude, but also a different shape. As a 
consequence, the wavelength- and position-dependence of the brightness distribution cannot 
be separated B{\,p) ^ /i(A)/2(p) and equation (||) truly contains a chromatic dependence 
of the amplification. This fact does not contradict the achromaticity of light deflection 
predicted by the general theory of relativity, it merely reflects the different appearance of 
the source star at different wavelengths. 

For most microlensing events the point-source limit in equation (^ gives a satisfactory 
description of the light curve. Extended source effects become important only in events with 
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a sufficiently low impact parameter (defined as the projected distance of closest approach 
of the lens and source on the sky). Expanding the amplification (^ in powers of the inverse 
of the lens-source separation a^^ irrespective of the Einstein radius of the lens (keeping the 
ratio e/cTo fixed) we obtain 



A{X,ao) 



aoJa^ + Ae^ 



1 







(6) 



Strictly speaking, this ctq ^ 1 expansion is convergent for ctq > 1 + V^- The first extended 
source correction term in the brackets is a product of two factors. The first one depends 
purely on the lens parameters (e, o"o), while the second one depends only on the source 
properties (B). The lens factor is a monotonously decreasing function of the lens distance, 
dropping from a value of Cq'^ /A when the source lies well within the Einstein ring (ctq <C e) 
to a value of Se^^a^^ when the source is well beyond the Einstein ring (ctq ^ e). The 
magnitude of the chromatic source factor is less than one, typically on the order of 0.4. 
Using this value we find that any lens with e ^ 4 will produce a 1% effect at ctq = 3. 
Therefore, in order to achieve a larger extended source effect, even a lens with high e 
has to approach the source closer than 3 source radii, beyond the region where the above 
expansion is validQ. An observed example of such an event is M95-30 (Alcock et al. 1997d), 
with a fitted value of the impact parameter p = 0.715 ± 0.003 and e = 13.23 ± 0.02. 

Note that we estimate the extended source effect in comparison to a point-source light 
curve with the same event parameters. Studying detectability requires comparison with 
a best-fit point-source light curve, which further lowers the upper limit for the closest 
approach. Gould & Welch (1996) estimate that using color observations the finite-source 
effect is detectable out to 2 source radii. 

We study the case of close approaches or source transits directly from equations (^) 
and (^. For evaluation of the involved surface integral we follow the method described 
in Paper I. The a~^ divergence of the amplification factor Ao{cr) at the point directly 
behind the lens can be avoided by integrating in lens-centered polar coordinates. The 
product of v4o(a") with the area element dT, = a da d(j) then leads to a finite integrand. For 
B{X,p) we use brightness profile data generated for each wavelength directly from model 
atmosphere calculations (see §2). We interpolate the computed center-to-limb data points 
by simple quadratic segments a + hp^ . The computational advantage of this approach is 



^ Lenses with e <^ 0.4 are too weak to trigger a microlensing alert in the current surveys. Lenses with 
0.4 < e < 4 require an approach closer than 3 source radii for a 1% effect. 
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Fig. 3. — Lower panel: light curves for two selected wavelengths (451 nm - solid line, 759 nm - 
dashed line) of a T=3750 K, log (7=0.5 model atmosphere (spectral resolution as in Figure §), 
lensed by an e = 13.23 lens with zero impact parameter. Upper panel: corresponding 
normalized surface brightness profiles. For the purpose of alignment with the lower panel 
the profiles are reflected to negative values of p. 

that the integral can then be analytically reduced to one dimension (for details see Paper 
I). Additionally, as the lensed flux is linear in the interpolation coefficients, it is sufficient to 
compute the light curve for a single wavelength and then by substituting coefficients obtain 
the other light curves for the same event. With this technique it is feasible to perform the 
computation for a large number of wavelengths and thus study spectral changes associated 
with microlensing. 



3.2. Light Curve Shapes and Low Resolution Spectra 
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Fig. 4. — Upper panel: low resolution unlensed spectrum of the T=3750 K, log5'=0.5 model 
atmosphere (flux in units of W m'"^ Hz^^). Lower three panels: amplification by an e = 
13.23 lens as a function of wavelength for lens positions ctq = 0, 0.72 and 0.96. Vertical 
dotted lines mark two regions with TiO bands affected by microlensing in opposite ways. 
These regions are enlarged with high spectral resolution in Figure & 



The effect of the brightness profile shape on the light curve is demonstrated in Figure ^ 
For both selected wavelengths the profile is normalized to unit total flux, in order to match 
the form in which it appears in the amplification formula (§). The light curves plotted 
in the lower panel have a zero impact parameter, thus the time t measured from closest 
approach in units of source radius crossing time is equal to the lens distance (| 1 1 = (Tq). For 
illustration we use the M95-30 lens with e = 13.23. From the figure it can be seen that the 
light curve traces the brightness profile shape - the flatter profile produces a flatter light 
curve. In addition, the light curve with lower amplification at the source center typically 
crosses the other at ctq ~ 0.7 and has a higher amplification at the limb of the source. We 
note that flatter profiles usually indicate formation higher in the stellar atmosphere. 
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Figure ^ illustrates the wavelength dependence of the amplification on a low resolution 
(R=500) spectrum from 401 to 775 nm. These results are obtained using the same model 
atmosphere as in Figure ^, with T=3750 K and logg = 0.5, lensed by an e = 13.23 lens. 
The model spectrum in the absence of the lens is shown in the upper panel. The lower 
three amplification plots correspond to lens positions ctq = (source center), 0.72 (closest 
approach of M95-30) and 0.96 (close to the limb). Comparison of the plots at the center 
and the limb further illustrates the relative change in amplification at different wavelengths 
seen in Figure ^ — the amplification curve at the limb is virtually a vertical flip of the 
curve at the center. The overall slope of the curves indicates that for ctq = the spectrum 
will appear bluer and for ctq = 0.96 redder than in the absence of a lens (see also Gould & 
Welch 1996; Valls-Gabaud 1998). 

The finer structure of the curves provides further interesting results. Sections of 
the curves with significant spectral variations of amplification indicate spectral regions 
particularly sensitive to microlensing. Observations of such regions during a transit event 
can then be used to test the depth dependence of the atmosphere model. While the 
resolution of this particular spectrum is too low to resolve individual spectral hues, broader 
molecular bands are clearly visible. An example of a sensitive feature is the TiO band 
system at 710 nm {A — X, 7-bands), as hinted by the bump in the amplification curve. 
Recent models by Jacob et al. (1999) also point to the extreme limb-darkening effects in 
this band system. It should be noted that not all bands are affected in the same manner — 
the C — X TiO system at 516 nm (a-bands) shows an opposite variation in amplification. 
The nature of the change in these bands is discussed in more detail further below. 

The three amplification curves are each plotted for convenience with a different vertical 
scale. For example, the relative amplification variation at cxq = 0.72 is much smaller than 
at either the center or the limb. A simple way to quantify the degree of this variation is to 
take the maximum, minimum and average amplification over the full range of wavelengths 
and plot the ratio {Amax — Amin) / Aaver as a function of lens position ctq- This ratio in 
fact serves as a measure of chromaticity - it drops to zero in the achromatic point-source 
limit. Results for three different atmosphere models (wavelength range 401-775 nm) are 
presented in Figure ^. While there is a difference in overall amplitude, the general character 
of all three curves is similar. As the lens approaches the source, the degree of chromaticity 
increases, reaching a local peak of ~ 3-5% just after crossing the limb. If the lens passes 
closer to the center, the chromaticity drops to < 0.5% at ctq ~ 0.7 (approximately the 
M95-30 impact parameter, by coincidence). This indicates that the light curves for all 
the studied wavelengths intersect approximately at the same point, as seen in Figure ^. 
A yet closer approach causes the chromaticity to increase to a peak value of ~ 8-11% 
when the lens is positioned at the source center. To summarize, the strongest spectral 
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Fig. 5. — Variation of amplification within the low resolution spectrum as a function of lens 
position for three model atmospheres: T=3750 K, log (7=0.5 (sohd line), T=3500 K, \ogg=l 
(dashed line) and T=4000 K, log (7=0.5 (dotted line). The Einstein radius of the lens is the 
same as in previous figures (e = 13.23). 

effects can be expected when the lens is aligned with the source center, with weaker effects 
near the limb. In between these two positions there is a drop in chromaticity nearly to 
zero (achromaticity). Although the limb effects are weaker, in principle they should be 
observable in any transit event. 

The results in Figure |^ were computed using the M95-30 lens with e = 13.23. The 
dependence of the results on the Einstein radius is explored in Figure ^ for the model 
atmosphere with T=3750 K and log (7 = 0.5. While the chromaticity predictably increases 
with e, its ao-dependence rapidly reaches a limit curve. In other words, the degree of lensing 
chromaticity cannot increase beyond this curve for a given atmosphere model and fixed 
spectral resolution. This behavior can be understood by studying the form of equation (^ 
during a close approach (as compared to the Einstein radius). In the limitQ e ^ ctq the 



^ Convergence requires also e > {ao + l)/2 . 
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Fig. 6. — Variation of amplification within the low resolution spectrum (401-775 nm) as a 
function of lens position for the T=3750 K, log (7=0.5 model atmosphere. Curves correspond 
to lenses with different Einstein radii in units of the source radius (from lowest curve): e = 
0.1, 0.2, 0.5, 1, 2, 5, 100. 
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is directly proportional to the Einstein radius e. The degree of chromaticity depends on 
amplification ratios and is therefore independent of e in this limit, as seen in Figure ^. We 
conclude that as far as spectral effects are concerned, the "only" advantage of microlensing 
by stronger lenses (e > 5) is an overall increase in flux. 

Figure § illustrates the capability of microlensing to resolve the depth structure of 
a given stellar atmosphere. Can the microlensing signature also be used to distinguish 
between different atmospheres? Figure |^ contains the results of Figure ^ together with those 
computed for the two other atmosphere models used in Figure ^; a cooler one (T=3500 K, 
\ogg=l) and a hotter one (T=4000 K, log5'=0.5). As expected, the results indeed vary 
substantially between the different models. In the 3500 K atmosphere, for example, the 
TiO a-bands at 516 nm show a stronger ampliflcation variation, and more importantly. 
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Fig. 7. — Microlensing by an e = 13.23 lens of three different model atmospheres: T=3500 K, 
log (7=1.0 (left column); T=3750 K, log5f=0.5 (central column) and T=4000 K, log(?=0.5 
(right column). Upper row: low resolution unlensed spectra (flux in units of W m^'^ Hz^^); 
central row: amplification for lens position ctq = 0; lower row: amplification for lens position 
o"o = 0.96. Vertical dotted lines mark the TiO a-bands (516 nm) and 7-bands (710 nm). 

the TiO 7-bands at 710 nm show the opposite variation than in the 3750 K case. The 
amplification variation in the 4000 K model is much weaker overall, because the TiO bands 
are less prominent at this higher temperature. The only prominent amplification feature at 
420 nm is in a spectral region with very low flux, hence it would be difficult to detect any 
spectral changes there. 



3.3. High Resolution Spectra and Line Profiles 
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Fig. 8. — Upper row: high resolution unlensed spectrum of the T=3750 K, log (7=0.5 model 
atmosphere (flux in units of W m~^ Hz~^). Lower two rows: amplification by an e = 13.23 
lens as a function of wavelength for lens positions uo = and 0.96. Left panels: spectral 
region near the 516 nm TiO a-bands. Right panels: spectral region near the 710 nm TiO 
7-bands. 



The results presented in Figures ^, |^, ^ and |^ illustrate the overall microlensing 
effects on low resolution spectra. However, studying the behavior of finer features, such 
as individual molecular bands or spectral lines, can provide stronger constraints on model 
atmospheres. First we return to the TiO bands mentioned above. In this particular case, 
even the low resolution results illustrated in Figure ^ indicate the different behavior of the 
two band systems. When the lens is at (Tq = 0, the 7-band system at 710 nm has a higher 
amplification in its region of highest absorption. As a result, the band system will appear 
weaker (shallower) than in the absence of the lens. With the lens at the limb, the highest 
absorption region has a lower amplification than the adjacent spectral region. The band 
system will therefore be more prominent than if the lens were absent. In contrast, the 
a-band system at 516 nm shows exactly the opposite behavior. It will be more prominent 
with the lens at the source center, and less prominent when the lens is at the limb. The 
different behavior of the 7-band system is due to the extension of its region of formation 
higher in the atmosphere. 

These effects are verified by the results of high spectral resolution (R=500,000) 
calculations shown in Figure ||. The first feature to notice in these plots, apart from the 
confirmed band system structure, is the high "fuzziness" of the amplification curves. This 
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Fig. 9. — Variation of amplification as a function of lens position for the T=3750 K, log (7=0.5 
model atmosphere lensed by an e = 13.23 lens. Low resolution spectrum 401-775 nm (solid 
line), high resolution 500-531 nm (dotted line) and high resolution 700-740 nm (dashed line). 

property indicates that the amplification varies highly within individual spectral lines. In 
fact the less "fuzzy" 710 nm region also has a lower density of spectral lines. Studying 
the amplification variation of both regions in Figure P, it comes as no surprise that while 
the shape of the dependence is similar, the degree of chromaticity is higher than in the 
full extent of the low resolution spectrum. This higher variation is due to the fine spectral 
structure, which is averaged out in the low resolution results. The degree of chromaticity 
for the 516 nm and 710 nm regions in high resolution again peaks at the source center (20% 
and 12%, respectively) and has a secondary peak at the limb (9.5% and 5.5%, respectively). 
For comparison, the center and limb values for the low resolution spectrum are 9% and 
4%, respectively. Note that the high values for the 516 nm region are partly caused by the 
single atomic spectral line at 526.95 nm Q 

A closer look at the fine structure in both sets of high resolution results (apart from 
the large-scale TiO band effect) suggests a significant correlation between the unlensed 



•^ Unfortunately this line is highly saturated at its core, where the amplification undergoes the high 
variation - the effect will therefore be difficult to detect because of the low photon flux at the line center. 
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Fig. 10. — Detail of Figure |^ resolving individual spectral lines. The spectral region shown 
is between 738.8-739.2 nm. 



spectrum and the amplification curve for o"o = 0. A blow-up of a 0.4 nm interval in 



Figure [T^ confirms the correlation - most individual lines have their cores amplified less 
than their wings at ctq = 0, while at aQ = 0.96 the cores are amplified more than the wings. 
Hence most lines become more prominent with the lens at the source center, and less 
prominent with the lens close to the limb than in the absence of the lens (ctq — >■ oo). There 
are fewer examples of lines with the opposite behavior (e.g. at 739.044 nm in Figure p^ ), 
with virtually no change (fiat amplification curve) or with more complex behavior. 

The characteristic correlation can be explained for typical absorption lines in the 
following way. Such lines form at a certain depth in the atmosphere, and gradually become 
less prominent away from the center towards the limb, as there is relatively less absorption 
at the limb. The center-to-limb variation curves at the cores of these lines are therefore 
flatter than in the adjacent continuum. The correlation is now directly explained by the 
discussion of light curve shapes above - when the lens is at the source center, the line 
cores are amplified less than the continuum, while with the lens at the limb the cores are 
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Fig. 11. — Equivalent width of Ha (left panel) and B.j3 (right panel) as a function of lens 
position, for the T=3750 K, log (7=0.5 model atmosphere. Individual curves correspond to 
different lenses (from lowest at o"o = 0): e = 0.1, 0.2, 0.5, 1, 2, 5, 100. 



amplified more than the continuum. 

The less common anti-correlation (or lack of correlation) of the unlensed spectrum and 
the amplification at source center can be expected for example in some emission features, 
high excitation lines, or in spectral lines affected by more complex physical effects, such as 
resonant scattering (Loeb & Sasselov 1995), overlapping lines, lines within molecular bands, 
etc. 

A straightforward method for studying the change of a particular line is to compute its 
equivalent width, which is generally given by 



W{ao) = A 



1 



JF{X,ao)dX 
JF,{X,ao)dX^ 



where Fc(A, ctq) is the continuum flux and the integrals are performed over a wavelength 
interval of width A around the line. As an example we present here results for the hydrogen 
Balmer lines Ha and H/5 in the same T=3750 K, log(7=0.5 model. The center-to-limb 
variation of their unlensed profiles is shown in Figure |l[ The dependence of their equivalent 



widths on lens position is plotted in Figure |TT| for a set of Einstein radii e. First we note 
that the change with e has the same character as in Figure |^. Any sufficiently strong 
lens (e > 5) will therefore affect the equivalent width in the same way. This result can 
again be traced to the fact that the equivalent width in equation (pi) depends on a ratio of 
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amplifications. This ratio is independent of e in the e 3> ctq hniit, as seen from equation (J^). 
During a microlensing event both Balmer hues behave in the typical manner described 
above. An approach of a sufficiently strong lens first causes the width of both lines to drop 
(by ~7%), as it gives higher weight to the limb where the lines are weak. Similarly, a lens 
positioned at the center, where both lines are strong, would increase their width (~19%) 
above the asymptotic value. Both the Ha and H/3 curves are similar, although those of H/? 
are more centrally peaked and have a broader minimum near the limb. 

We compared our results to those of Valls-Gabaud (1998) for H/3 with the same event 
parameters^. The general character of the change with lens position is the same - decrease 
of the equivalent width (EW) at the limb and a higher increase closer to the source center. 
However, Valls-Gabaud's results show a two to three times lower relative change of the 
EW, the actual factor depending on the lens position. As a result, the predicted ratio of 
EW increase (close to the center) to peak EW decrease (at the limb) is also lower than our 
computations show. Furthermore, the achromaticity point appears to be closer to the limb 
(around ctq = 0.8) than in our analysis. As mentioned in §2, assumptions such as linear 
limb darkening degrade Valls-Gabaud's models too much to give realistic predictions of 
equivalent widths, particularly in the case of red giants. This comparison further illustrates 
the sensitivity of microlensing to the atmosphere structure of the resolved source. 

A detailed study of high resolution spectral changes in different model atmospheres and 
comparison with the M95-30 observational results will be presented elsewhere (Heyrovsky 
& Sasselov 2000b). 



4. Conclusions 

Microlensing events in which a point-mass lens with Einstein radius e > 4 (in source 
radius units) approaches the source star within three source radii already exhibit ~1% 
deviations from point-source light curves. The deviations further increase for closer 
approaches. These finite source effects depend on the surface brightness distribution of 
the source star. Observations of such effects, particularly during source transit events, 
can therefore be used to resolve and study the structure of the stellar disk. As the 
brightness distribution depends on wavelength, the usual assumption of microlensing 



"^ We neglect here the earher resuhs m Vahs-Gabaud (1996), which are mconsistcnt with Valls- 
Gabaud (1998). 
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achromaticity breaks down in these cases. Spectral observations can then be used to probe 
the depth structure of the atmosphere of the source star. The most promising sources 
for such studies are red giants in the Galactic bulge, among other reasons due to their 
large size, intrinsic brightness and fairly high abundance. These circumstances are in fact 
particularly fortuitous, as the structure of cool red giant atmospheres is currently very 
poorly observationally constrained. 

The degree of microlensing chromaticity can be measured by the relative spectral 
variation of the microlensing amplification. Thus defined, chromaticity as a function of lens 
position is found to have two peaks - primary peak at the source center and secondary peak 
close to the source limb, with a dip to achromaticity in between. In the usually observed 
type of events, in which the Einstein radius e > 5, chromaticity is independent of e. Its peak 
values for a low resolution optical spectrum are typically 10% at the center and 4% at the 
limb (both values increase with higher spectral resolution). Any such transit event should 
therefore exhibit at least the 4% effect, which should be readily observable under favorable 
conditions. In future potential source- crossing events it is therefore important not to miss 
the time when the lens is close to the limb, as was unfortunately the case in M95-30. 

During a red giant transit event the overall spectrum appears redder when the lens is 
at the source limb, and bluer if the lens comes close to the source center. Individual spectral 
lines respond to microlensing in different ways, because of their different center-to-limb 
variation, which reflects their different depth of formation in the star's atmosphere. Most 
absorption lines will typically turn weaker (lower equivalent width) when the lens is close 
to the limb and more prominent when the lens is near the center. Some emission lines, high 
excitation lines, overlapping lines, lines within molecular bands and lines affected by more 
complex physical effects can exhibit opposite or generally different behavior. As a result, 
broad molecular bands visible even on a low resolution spectrum of a given source may 
behave in different ways. The pattern of their changes depends on physical parameters, 
composition and structure of the atmosphere. Observations of these bands can therefore 
be used to constrain the global parameters, and to check theoretical assumptions used for 
constructing model stellar atmospheres. High resolution spectra resolving individual lines 
would, however, also enable studies of the depth structure of the star's atmosphere in detail. 
Particularly in the case of red giants, gravitational microlensing therefore provides a unique 
new tool for stellar physics. 

During M95-30, the only well-documented point-mass microlensing transit event up 
to now, eleven low resolution and three high resolution spectra were taken in addition 
to photometry in two colors (Alcock et al. 1997d). In a forthcoming paper (Heyrovsky 
& Sasselov 2000b), we will use the methods and findings described here together with 
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the observational data to study the lensed giant and construct an adequate model of its 
atmosphere. 

We thank R. Kurucz for his kind help with opacities and molecular data. We also 
thank an anonymous referee for comments that helped improve the manuscript. DS is 
grateful to I. Shapiro for his support during part of this study. 
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